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THE SEPARATED TURBULENT BOUNDARY 


LAYER OVER A WAVY WALL 
NASA Grant No. NSG 1208 

Final Report 
covering the period 
July 1, 1975 - September 30, 1977 

by 

A. Polak and M.J. Werle 
University of Cincinnati, Cincinnati, Ohio 

This research was concerned with the development of a 
prediction method for calculating detailed distributions of 
surface heating rates and pressures over a corrugated wall, 
due to flow of a two-dimensional supersonic turbulent boundary 
layer. Primary attention was given to the turbulent boundary 
layers which are thick relative to the wave ampiicude and at 
the same time the ratio of the wave amplitude to the wave 
length is sufficiently high so that a strong interaction develop 
between the viscous and inviscid regions. This manifests itself 
in the variation in surface pressure which may be an order of 
magnitude different from the one calculated by purely inviscid 
methods. Also, under these flow conditions, separation regions 
may develop between the peaks of the waves. This means that 
the nature of the problem is such that classical boundary layer 



methods are not suitable to effect the solution, simply be- 
cause the interaction is a first order effect. 

To solve this problem effectively, an interacting boundary 
layer approach was employed. Its essence is that the pressure 
distribution is not prescribed by the body surface profile, but 
calculated simultaneously with the viscous flow. A time-like 
relaxation method using an implicit finite-difference numerical 
code was adopted. The turbulence was represented by the two 
layer eddy viscosity model of Cebeci and Smith. An attempt was 
made to account for the surface curvature effect on the turbulent 
structure by employing a correction to the eddy viscosity model 
proposed by Bradshaw. The present predictions were compared 
with independent experimental and theoretical results. The pre- 
dictions, in terms of pressure and heating levels, compare reason- 
ably well with the measured data when the curvature correction is 
included . 

The complete details, concerning the formulation of the 

problem, the numerical method and the discussion of the results 

for test cases with Mach numbers of 2.5 and 3.5, Reynolds numbers 

Re = 10.8xl0^/mm and 32.4xl0^/mm, wall-to-stagnation temperature 
00 

ratio T^/T-j = 0.4 0 and 0.81, and ratio of wave amplitude to wave 
large from 1/96 to 1/12 are presented in a self-contained 
University of Cincinnati technical report. Report No. AFL 77-11-36 
(Reference 1). The results of this research effort were also presented 
in a paper at the ASME Winter Annual Meeting, Atlanta, Georgia, 

November 1977. This paper appears as an ASME publication 77-WA/HT 
-28 (Reference 2). While this work was in progress partial results 
of this research were presented at the 5th U.S. Air Force/Federal 
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Republic of Germany Data Exchange Agreement Meeting, Wright 
Patterson Air Force Base, Ohio, 28-29 April 1976 (Reference 3). 

The present investigation lead also to two side efforts, 
carried out under partial support of the NASA Grant No. NSG 1208. 
These efforts were aimed at achieving higher accuracy and effic- 
iency in the numerical schemes, applicable to problems of present 
interest. The first of these studies was concerned with developing 
a numerical scheme for solving viscous flow problems with regions 
of predominantly inviscid flow (like the outer region of very 
thick turbulent boundary layer). This effort resulted in a new 
second order accurate exponential, box scheme. It was shown that 
this method provides a definite advantage over other numerical 
schemes. The details of this work were published as a University 
of Cincinnati technical report, Report No. AFL 77-9-34 (Reference 
4) and also resulted in a Master of Science Thesis. A Technical 
Note summarizing this work was submitted for publication to the 
AIAA Journal. A copy is attached herein as Appendix A. 

The second study aimed at achieving a higher efficiency in 
calculating laminar and turbulent boundary layers. The fourth 
order cubic spline method developed recently by Rubin was combined 
with the Davis Coupled Algorithm for solving simultaneously the 
boundary layer equations. The details of this work, with appli- 
cation to a set of model flow problems and to the incompressible 
laminar boundary layer equations will appear as a University of 
Cincinnati Technical Report (Reference 5) in December 1977. This 
work appears also as a Master of Science Thesis; the abstract of 
this Thesis is attached herein as Appendix B. 
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In conclusion, it car. be stated that the major objective of 


tne work under the present NASA Grant No. NSG 1208 was fulfilled. 

A method capable of handling strongly interacting flows was success- 
fully adapted to solve the problem of flow of thick turbulent 
boundary layers over a wavy wall, including multiple separation 
regions. A supportive effor + resuir.ed in developing numerical 
schemes, which in the future could enhance the accuracy and 
efficiency of the present state of the interacting numerical 
algorithm. 

The cumulative list of publications resulting during the 
course of investigation under the support of the NASA Grant No. 

NSG 1208 is attached under the heading List of References. 
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APPENDIX A 

A NUMERICAL METHOD TCP. SOUMDARU LAVERS WITH BLOWING 
- THE EXPONENTIAL SOX SCHEME* 


T.M. El Mistikawy and M.J. Aerie 

Department of Aerospace Engineering and Applied Mechanics 

University of Cincinnati 
Cincinnati, Cr.io 


INTRODUCTION 

Solutions of boundary layer flews in the presence of massive 
injection through a pore us surface are complicated by the multi- 
structure nature of the flow field. The flow can be divided 
into two distinguishable regions (see ?.ofs. 1-41? a basically 
in vise id turning region near the blowing surface that occupies 
most of the boundary layer and a thin viscous region in which 
diffusion effects olav a sicnificar.t role in adjusting the 

i 

boundary layer flew to the free stream conditions. Similar 
flew structure and complication occur in a wide class of 
viscous interaction flow problems f?efs. 5 through where a 
thin viscous layer actively induces a significant displacement 
effect to a relatively thick mviscid flew layer. The change 
ir. the character of the flew from the vis ecus to in viscid 
regions is difficult to accurately capture numerically with a 
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single finite difference representation giving rise to the 
introduction of various schemes to eliminate numerical ‘ • 
"instabilities" ranging from variable grid spacings (?.ef. 10) 
to artificial viscosity concepts (?efs. 11 and 12). 

The present study takes a different approach and presents 

a new numerical scheme that should be useful for solution of a 

wide class of inviscid/vi3cous flews. The approach is based on 

two points; first, that the difficulties encountered \ previous 

schemes can be traced to a cusr epresentaticn of convection 

effects normal to the principle boundary layer direction (see 

?.ef. 13) and, second, that exponential difference operator 

concepts (P.efs. 14 through 13) can be combined with Keller's 

box scheme approach (Kef. 19) to produce a stable second order 

accurate finite difference scheme for such ccnvecticn-dif fusion 

•• 

problems. The technique developed here is demonstrated 
through application to the solution of the Falkr.er-Skar. ecuatior.3 
for self similar boundary layer flew with massive blowing at 
the surface. 


gcvz?:;t::g z:tat:c::s ant scivt:::; 


The governing equations are taken tc be the Falknsr-Skar. 
equations with blowing and giver, here in terms of the Gcrtlar 
variable n (see Per. 20) as 


* F * 0 





(la) 

(lb) 




with 


F (0 ) - 0 ; V{0) - V ; (2a) 

w 

and F - 1 as n - - * * (2b) 

where V is related to the physical injection velocity, F represents 
the normalized longitudinal velocity component, and 3 the inviscid 
pressure gradient parameter. 

As in nearly ail numerical solution methods, the nonlinear 
nature of equation (lb) is addressed indirectly. An iterative 
approach is used to first write this in the form 

F - (a+b) F + lb F ■ c (3) 

d n n 

where a, b are obtained frcm a previous iteration and are given 


as 

a 

5 V/2 + 

/(V/2 ) 2 + SF 

(4'a) 


b 

» V/2 - 

/(V/2) 2 r 3F 

(4b) 

with 

c 

- - 3 


(4c) 


Whereas classical finite difference techniques see!*: solutions 
of equation (3) using Taylcr series to generate difference 
approximations to derivatives over a small grid distance, a 
different approach is employed here. Following the 
lead of references 14-13, the coefficie 

a ra >- a era 1 1 ' A /-3 s 


5 Of 

equation (2) are 

, an ce 

and the resulting 

nee 


the 

O w « •* « * 

os 1 

and 2 ia given as 




!T rgr&rx (\\ Iff 








a.n b.n 

* 1,2 A 1 * S 1 e -1 


( 5 ) 


where f. ■ c,/a.b, and the ccef f icients a. and b, are evaluated 


X x x 


at the rr.idpcir.t3 of the grid indicated. Use cf equation (5) 
along with a similar expression for F, 2 ia the interval between 
poi-ts 2 and 3 and the two continuity conditions that F ?> ■ F-._ 


ana 


'2t 


2 - 


maxes it possible to eliminate the arbitrary 


constants A, and 3. . The resulting difference equation is 

a l F ! " (S 1 + 3 2 )F 2 + a 2 F 3 * (a l “ 3 l )f l + (a 2 “ i 2 )i 2 


where 


i. * (a.-L.)e 


-a-in. 


b-j i n 1 - to'i* 

i /Y i ; »2 ” (a 2 -b 2 )e “ Vy 2 


[a : - (1-Y 1 )b 1 ]/Y 1 ; 3 2 - (a : (1 -y 2 ) -fc-J/Y 

- (a,-b, ) in. 


1 . 


y 2 ■ 1 - e 


“ (a 2 -b ? ) i " 2 


( 6 ) 

(7a) 

(7b) 

(7c) 


Equation (6) i3 a three point difference approximation to 
equation (3) and, as such, produces a tridi agonal set of algebra! ■ 
equations that are straightforward to solve using standard 
techniques. In addition, this difference relation possesses 
all the favorable properties of the exponential schemes and 
the "2-poir.t" or "box-scheme” approaches. Equation (5) is 
found to be always diagonally dominant for ab < 3, it is second 
order accurate in in. throughout the diffusion region, and most 


* The use of two points to generate the basic relations for this 
scheme aligns it with the Bex Scheme techniques of Feller. 

Pef. (19); hence the term "exponential box-scheme” emerges. 
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importantly, it automatically and snocthly switches to a second 
order accurate windward difference scheme as the normal convection 
velocity, V (as represented by the term a+b) , becomes large 
negative cr positive. 

The companion numerical solution of the continuity equation 
is obtained substituting equation (5) directly into equation (la) 
and integrating to obtain the local relation 

a. An, b,An, 

v 2 - v i + A^l-e 1 1 )/a 1 + B^l-e X A )/b 1 - (8) 

where A. and 3^ are obtained in terns of f . and using equation (5) . 

PES'JLTS AMD SISCVSSICM 


The exponential lex (Z3) scheme described above hts been 
applied ahd verified for a wide range of test cases (see Pef. 13) 
with typical results presented here. Typical longitudinal 
velocity profiles are shewn in Figure 2 for blowing at an axi- 
sytanetric stagnation point (3 • O.j). The massive blowing case 


(V * 25) is similar to that presented by Liu and Chiu ( P.e f . ID) 
w 


except that in that study, compressible heating effects were 
accounted for due to a cold wall boundary condition. The present 
case corresponds to the case of an adiabatic wall temperature 
and was chosen to allow isolated study of the Z3 scheme for 
the momentum and continuity equations alone. 

In an attempt to demonstrate the multistructure nature of 
the flow field. Figure 2 shews the individual component con- 
tributions to the governing momentum balance equation ,1b) in 




;r* tic** renen 


• i /■* » 


(4! 


< 55). Clearly the 
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flow represents an inviscid process beicw - < 43 where the 
normal convection and pressure gradient effects totally balance 
one another. Only above this region does the diffusion effect 
play i significant ro-e. With so much of the flow basically 
inviscid in nature, it should not be surprising to find that 
diffusion- type numerical schemes (such as the classical three 
point central difference schema) have difficulty representing 
this flow. Such difficulties usually manifest themselves in 
terms of non -mono tonic (jagged) velocity distributions (cr local 
"wiggles'*). Thi3 point is discussed further in Ref. 12 with the 
conclusion that this difficulty is largely due to a failure of 
central difference schemes to adequately medal the normal 

convection effects. It is shown in Pef. 13 that a "windward" 

/ 

representation of the convection terms is necessary to properly 
represent such flows. The exponential box scheme developed here 
automatically and smoothly switches from a second order accurate 
diffusion- type scheme in diffusion regions (n > 43) to a 
second order accurate windward scheme in normal convection dominate 
regions (n < 43), thus encountering no numerical difficulties cr 
"instabilities" even for a constant grid mesh as per ?.<= f . 10. The 
accuracy and reliability of the current approach is attested to 

in Figure 4 where the results cf a typical step size study i3 shewn 

*• , 

for 3 points in the critical region cf the flow-field, 45 ^ n 
Here the E2 scheme is seen tc maintain second order accuracy ' i.e. 
produce a straignoline variation with u- “) for a uniform step 
size An as large as 2 while a central difference scheme (CD) 
rapidly loses its second order nature. Reference 13 presents a 
more detailed comparison and srtev/s the current approach tc be 
superior tc four ether difference methods as well. 
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APPENDIX B 


AN APPLICATION OF FOURTH ORDER SPLINE COLLOCATION 
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ABSTRACT 


A study and application of the fourth-order eplino 
collocation procedure, numerical solution of boundary laver 
like differential equations, is presented. A simple 
inversion algorithm (represting a fourth-order version of 
the Davis Coupled Scheme) for the simultaneous solution 
of the resulting difference equations is given. Particular 
attention is focused on the boundarv condition representation 
for the spline second-derivative approximations. Solutions 
using the spline procedure, as well as the three-point 
finite difference method, are presented for several model 
problems in order to assess and improve the spline numerical 
scheme. Application of the resulting algorithm to the 
incompressible laminar self similar boundary layer 
equations is presented. 


